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Abstract 

The present status of the renormalization-group (RG) predictions 
on the superfluid density p s and the specific heat C + and C~ near the 
A- line of 4 He is briefly reviewed. Particular attention is given to uni- 
versal amplitude ratios related to these quantities. The goals of a new 
theory project are presented that involves higher-order calculations of 
the amplitude functions of p s ,C + and C~ which are of fundamental 
importance for a test of the universality predictions of the RG theory. 

1 Introduction 

One of the fundamental achievements of the renormalization-group (RG) 
theory of critical phenomena is the identification of universality classes in 
terms of the dimensionality d of the system and the number n of components 
of the order parameter M. Specifically, RG theory predicts that, within 
a given universality class, the critical exponents, certain amplitude ratios 
and scaling functions are universal quantities that do not depend, e.g., on 
the strength of the interaction or on thermodynamic variables (such as the 
pressure). The superfluid transition of 4 He belongs to the d — 3, n = 2 uni- 
versality class and provides a unique opportunity for an experimental test 
of the universality prediction by means of measurements of the critical be- 
havior at various pressures P along the A-line T\(P) (and along the A-line of 
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3 He — 4 He mixtures). Early tests have been performed by Ahlers and col- 
laborators and consistency with the universality prediction was found within 
the experimental resolution ||. At a significantly higher level of accuracy, 
the superfluid density p s and the specific heat C + and C~ (or, equivalently, 
thermal expansion coefficient (3^) above and below T\(P) are planned to be 
measured in the Superfluid Universality Experiment (SUE) || under micro- 
gravity conditions or at reduced gravity in the low-gravity simulator Q. As 
demonstrated recently ||, this would allow to perform measurements up to 
\t\ ~ 1CT 9 in the reduced temperature t = (T — T\(P))/T\(P). 

Important steps towards an improved universality test can already be per- 
formed, even in the presence of gravity, by new ground-based measurements 
provided that non-universal and universal effects are properly separated in 
a nonlinear RG analysis of the data [FJ. To extract the leading critical ex- 
ponents of p s and C ± from the experimental data and to demonstrate their 
universality at a highly quantitative level requires detailed knowledge on 
certain universal amplitude ratios and on non-asymptotic corrections to the 
asymptotic power laws. These quantities can be calculated by means of field- 
theoretic RG methods applied to the y? 4 model fl], §]. 

In the following we point to a serious lack of quantitative knowledge in the 
theoretical literature 0] on these amplitude ratios and correction terms. We 
also discuss the implications of very recent two-loop results for p s [|7| and 
shall summarize the goals of a new theory project || on higher-order RG 
calculations of the amplitude functions of p s and C^. 

2 Universal amplitude ratios 

In the previous analyses of experimental data for p s and C (or (3^) the 
following representations were employed 0, [|, |9| 

p s = k B T x (m/h) 2 A p (l + h |*| )\t\^ (1 + a p |t| A ) , (1) 
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C + = — \t\~ a (l + at\t\ A + E + \t\) + B , t>0 (2) 



a 



CT = — \t\- a (l + a'\t\ A + E~\t\) + B , t<0 (3) 



where m is the helium mass. These representations contain 10 nonuniver- 
sal amplitudes A p ,ki,a p ,A ± ,af,E ± ,B which need to be adjusted at each 
pressure. Their values vary by 15 — 30%, typically, along the A-line. The 
exponents a and ( = v{d — 2) = (2 — a)(d — 2)/d and the Wegner ex- 
ponent A are predicted to be universal. The existing theoretical predic- 
tions on these exponents |]RJ are based on field-theoretic calculations to five- 
loop order [11| and Borel resummation. The recent experimental result || 



a = —0.01285 ± 0.00038 is consistent with but more accurate than the RG 
estimate In addition, RG theory predicts the amplitude ratios 



A+/A-, (A~) 1/d /A p , a-/a p , a+/aj (4) 

to be universal Knowing their universal values for the d = 3, n = 2 
universality class would impose significant constraints on the analysis of the 
experimental data and would thereby improve the reliability and precision of 
the experimental results regarding the universality of a and (. 



The present theoretical knowledge of the ratios (4), however, is based only on 
low-order (mainly 1- and 2-loop) field-theoretic calculations |]J which imply 
an uncertainty at the level of at least 10 — 30%. For example, even the sign of 
the ratio a~ / a p was incorrectly predicted by the two-loop e(= 4—d) expansion 
[0. This issue was resolved by a (low-order) calculation within the d = 3 
field theory |]I3" |. The present status of the experimental and theoretical 



values of the universal ratios (4) is summarized in Table 1 (without error 
estimates). 
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Table 1 Universal amplitude ratios for the n = 2, d = 3 universality class 











(A-)V3 Mp 


a c /a p 


a>t/ a c 


Experiment 


1.067 


0.85 


-0.068 


1.03 


Experiment flUfl 


1.088 






0.85 


d = 3 Field Theory @ 


1.05 


0.78 


-0.045 


1.6 


e— Expansion 


1.029 


1.0 


1/6 


1.17 



As noted in Ref. 13, there is a large uncertainty of the theoretical value 
for <ic/a~. The d = 3 field-theoretical values for (A - ) 1 / 3 / A p and a~/a p are 
based on one- loop results [O]. Very recent d = 3 two- loop results M, to be 



discussed in the subsequent section, imply considerably different values [15 



thus indicating a considerable lack of reliability of low-order results for the 
amplitude ratios (4). 
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3 Recent results on p s in two-loop order 



In an effort to improve the theoretical prediction on these amplitude ratios, 
a two-loop calculation of p s was recently carried out by Burnett, Strosser 
and Dohm within the (p A field theory in d = 3 dimensions 0. The superfluid 
density can be calculated via Josephson's definition (16 



d 

p s = k B T(m/ri) 2 \i)o\ 2 —XT{k)' l / k= Q 



(5) 



where Xt(^) is the transverse susceptibility at finite wave number k and ipo 
is the order parameter (Bose condensate wave function). Within the d = 3 



<y2 4 field theory the representation (5) becomes |7|, [13 



(6) 



where 

= &\t\-»(l + aj\t\ A + ---) (7) 



is an appropriately defined correlation length. The amplitude functions fy(u) 
and fr{u) depend on the effective renormalized (p 4 coupling u = u^Z 1 ) which 
for — > oo approaches the fixed point value u(0) = u* = 0.0362 [|l^]. These 
functions are plotted in Fig. 1 in zero-, one-, and two-loop order. Their fixed 
point values fy(u*) and fr(u*) determine the asymptotic amplitude A p in 
(1) while their derivatives at u = u* contribute to the subleading amplitude 
a p . 
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Fig. 1: Amplitude functions ufy(u) and fr(u) for the order parameter and 
for the k 2 part of the inverse of the transverse susceptibility in zero-, one- 
and two- loop order vs the renormalized coupling u. The curves terminate at 
the fixed point u* = 0.0362. /From Ref. 7. 
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While the two-loop correction to fx is remarkably small (about 1%) at the 
fixed point u* = 0.0362, the two-loop contribution to fy(u) is about 10 %, 
thus indicating a considerable uncertainty of low-order perturbation theory. 
The uncertainty is even larger for the derivative of f^(u) at u* which implies 
a correspondingly large uncertainty of a p in (1) and of the ensuing amplitude 
ratio a~ /a p . A similar uncertainty exists with regard to the amplitude func- 
tion F_ (u) of the specific heat C~ below T A |13| determining A" and a~ 



in (3) as well as with regard to the amplitude function G(u) of the helicity 
modulus T which is proportional to p s . 



4 New theory project 

In order to significantly reduce the uncertainty of the theoretical predictions 
on the amplitude ratios (4) it has been proposed H] to perform new higher- 
order field-theoretic RG calculations and Borel resummations of various am- 
plitude functions. In general, the results of such calculations depend on the 
specific renormalization scheme employed. In view of various applications, 
the method of the minimally renormalized (p A field theory at fixed dimension 
d — 3 [171 [T3|, [T7| appears to be most appropriate. The following quantities 
need to be calculated: 

(i) The additive renormalization constant A(u, e) of the specific heat and 
the related RG function B{u). At present, A(u,e) and B(u) are known 
only in two-loop order |T7]]. For a three-loop calculation within a dif- 



ferent renormalization scheme see [20]. 



The amplitude function f^(u) of the square of the order parameter 
appearing in Eq. (6). For n = 1 it is known up to five- loop order and 
in Borel-resummed form |18|] but for n > 1 only up to two-loop order 



(iii) The amplitude function F_(u) of the specific heat below T\. For n = 1 
it is known up to five-loop order and in Borel-resummed form [|l8|] but 



for n > 1 only the two-loop result is known at present [13]. (Even for 
n = 1 the five-loop result needs to be revised on the basis of higher- 
order calculations for B(u) according to item (i) above.) 
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(iv) The amplitude function F+(u) of the specific heat above T\. Previously 
this function has been computed in Borel-resummed form for n = 1, 2, 3 
on the basis of six- loop results |0|. But this calculation needs to be 



revised on the basis of higher-order calculations for B{u) according to 
item (i) above. 

(v) The amplitude function f[(u) of an appropriately defined correlation 
length £_ below T\. This correlation length performs the task of absorb- 
ing logarithms in the four-point coupling uq in the bare perturbation 
series in d = 3 dimensions |7], 13J, |18|j . For n = 1, £_ can be taken to be 



the ordinary correlation length of three-dimensional Ising-like systems 
below T c . For a three-loop calculation of £_ for n = 1 within a different 
renormalization scheme see 



(vi) The amplitude function /t(w) [see Eq. (6)] of the k 2 part of the inverse 
of the transverse correlation function Xt(&)- At present this function 
is known in two- loop order [0]. Equivalently (and preferably), the am- 
plitude function G{u) = 47t/^(m)/t(m) of the helicity modulus T [|7| 
should be calculated in higher order. 

Higher-order calculations of G(u) or fx{u) would be more demanding than 
those of fip and F_ because they involve perturbation theory at finite wave 
number k. In view of the smallness of the two-loop correction to /t(w) (see 
Fig. 1) it has been argued [|7| that this low-order result for friu) can be 
considered as rather reliable, presumably within a few percent. In a first 
step of future calculations it will therefore be sufficient to confine the project 
to the quantities (i) - (v) noted above. 



Once the functions B(u), f^(u), F_(u), F + (u) and f^(u) are known at an ac- 
curacy of a few percent it will be possible to determine the amplitude ratios 
(4) at the same level of accuracy. It will be advantageous to go beyond the 
asymptotic representations (1) - (3) and, instead, include the entire Weg- 
ner series, i.e. higher-order terms |t| nA , n = 2,3,- • •, without introducing 
new nonuniversal parameters into the analysis. For the strategy of such a 
nonlinear RG analysis see Refs. 6 and 17. This may imply that the phe- 
nomeno logical analytic corrections ~ k\\t\ and ~ E \t\ in the representations 
(1) - (3) are unnecessary or will turn out to be much smaller than determined 
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previously [g §, g|. 



The success of this theory project depends on the order of perturbation theory 
up to which the calculations can be carried out. Preliminary considerations 
by Larin ]22| indicate that the amplitude functions /</,(«), F_{u) and /g~(u) 
can be determined up to four-loop order, and the RG functions A(u, e) and 
B{u), on the basis of previous work |nj, up to five- loop order. Such results 



will have an important impact on a future test of the fundamental RG pre- 
diction of critical-point universality. 
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Note added 

Very recently a five-loop calculation of A(u,e) and of B(u) has been per- 
formed for general n and a Borel resummation of B(u) has been carried out 
for n = 1,2,3 p3| j24|. The Borel resummed function B(u) differs from the 
two- loop result B{u) = ~ + 0(u 2 ) by less than 1 % at the fixed point it*. 
For n = 2 the result is B(u*) = 1.0053 ± 0.0022. Furthermore, a three-loop 
calculation of the amplitude functions f^{u) and F_{u) in three dimensions 
below T c has been carried out for general n > 1 []25]. [26| . 
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